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ABSTRACT 

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution 
method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter 
Telescope (BLAST). Compared to the direct Fourier transform method of deconvolution, the L-R 
operation restores images with better-controlled background noise and increases source detectability. 
Intermediate iterated images are useful for studying extended diffuse structures, while the later iter- 
ations truly enhance point sources to near the designed diffraction limit of the telescope. The L-R 
method of deconvolution is efficient in resolving compact sources in crowded regions while simultane- 
ously conserving their respective flux densities. We have analyzed its performance and convergence 
extensively through simulations and cross-correlations of the deconvolved images with available high- 
resolution maps. We present new science results from two BLAST surveys, in the Galactic regions 
K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. 

We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud con- 
taining K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, 
we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), 
bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diag- 
nostic tool to estimate the evolutionary stages of the clumps. There are close relationships between 
dust continuum emission and both 21-cm radio continuum and 12 CO molecular line emission. 

The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong 
spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 
12 K has been obtained for the central filament. We report physical properties of ten compact sources, 
including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact 
sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert 
mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless 
clumps has also been found in the IC 5146 region. 

Subject headings: submillimeter — stars: formation — ISM: clouds — balloons 
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1. INTRODUCTION 

Th e BLAST flight of 2 005 (BLAST05; IPascale et all 
120081: iChapin et all 120081 ) conducted both targeted and 
unbiased surveys in some of the important star forming 
regions of the Galactic Plane. These were conducted si- 
multaneously at the submillimeter wavelengths of 250, 
350, and 500 fim. BLAS T was a precursor to Herschel 
observations with SPIRE ([Griffin et al.ll2010h . One of the 
science goals of BLAST surveys was to study the early 
stages of massive star formation and identify their pre- 
cursors in order to describe a complete sequence of evolu- 
tion. In order to realize this ambition, a deep survey com- 
plemented by high spatial resolution is required. With a 
2-m telescope, BLAST05 optics were designed to produce 
dust emission maps at a diffraction-limited resolution of 
full width at half maximum (FWH M) 40, 58, and 75" a t 
250, 350, and 500 fim, respectively (|Pascale et al.ll2008l) . 
Unfortunately, due to some uncharacterized optical prob- 
lem, the BLAST05 Point Spread Function (PSF) was 
ano malous, producing images at a reso lution of about 
3^3 ([Truch et al.ll2008[|Rov et al.ll2010b| ). 

Most of the resulting maps contained clear imprints of 
the corrupted PSF (hexagonal geometry), indicating the 
presence of compact sources with angular sizes compa- 
rable to the diffraction limit. Measurements of the flux 
density of a few isolated bright sources can be carried out 
by aperture photometry or by fitting the PSF. However, 
because of the large beam area of the corrupted PSF, the 
sensitivity to compact point sources in the presence of 
cirrus noise suffers, and the detection of crowded sources 
is challenging. This leads to an underestimate of the 
source density and produces an incomplete source list. 

In order to improve the image resolution to near the 
diffraction limit, a direct Fourier transform method of 
deconv olution was perform ed in analyzing the Vulpecula 
region ([Chapin et al.ll200"8[) . In this method of deconvo- 
lution, compact sources in the processed map were en- 
hanced to 1.'3, 1.'6, and L'8 at 250, 350, and 500 fim, 
respectively. However, the deconvolved maps were con- 
taminated by ripples which propagated across the entire 
map. The generation of ripples is an intrinsic conse- 
quence of the direct Fourier inversion method when the 
PSF varies across the field and is not perfectly known. 
The amplitude of the ripples is particularly high near 
bright sources. Strong interference patterns are also ob- 
served in the vicinity of multiple bright sources. This 
makes fainter compact sources hard to detect; further- 
more, artificial spatial features are introduced in diffuse 
structures. 

To get the best science return from the BLAST maps, 
we consider implementation of the L-R method of de- 
convolution (|Richardsonlll972l : lLucylll974t ). The L-R de- 
convolution method is an iterative way of estimating the 
sky brightness distribution, conserving flux at every suc- 
cessive step. Fortunately, maps from the BLAST 2005 
flight have a high signal-to-noise ratio and are well sam- 
pled with 15" pixels, which is particularly important for 
the performance of L-R deconvolution. 

Implementation of L-R deconvolution on the BLAST 
maps, for example Cas A (iSibthorpe et al.|[2010f). Aquila 
(iRivera-Ingraham et al.ll2010f) . and Cvgnus X (Roy et al.l 
I2010al) . has remarkably enhanced the compactness of 
point sources to about 1' effective resolution. The other 



benefits of the L-R scheme are: 1) preservation of spa- 
tially diffuse structures; 2) suppression of background 
noise; and 3) resolving blended sources. As an example, 
in the Cygnus X region the L-R map led to a source cata- 
log containing 184 sources, significantly greater than the 
130 produced using the direct Fourier inversion method. 

A basic goal of this paper is to demonstrate the perfor- 
mance of the L-R algorithm in deconvolving BLAST05 
maps. We do this via simulations and using "ground 
truth," i.e., comparing the details in BLAST deconvolved 
images with available high-resolution maps, thereby es- 
tablishing the reliability of detected structures, both 
compact (§ 13.21) and more diffuse (§ 13 . 3[) . In those sub- 
sections we focus primarily on the operational conse- 
quences of deconvolution for identifying compact sources 
and defining diffuse structures. We subsequently present 
new science results made possible by this deconvolution, 
for the star-forming regions K3-50 (§ d]) and IC 5146 
(§0. We identify protostellar clumps and characterize 
their physical properties such as temperature (T), mass 
(M), and luminosity (L), by constraining spectral en- 
ergy distributions (SEDs) using multi-wavelength data. 
We have also investigated the properties in the "Northern 
Streamer," a cold molecular ridge in the IC 5146 field. 

2. OVERVIEW OF DECONVOLUTION 

Deconvolution is a common technique used by the as- 
tronomical community for enhancing image quality. In 
the imaging process, various factors are responsible for 
the distortion of the observed sky map. In the tradi- 
tional way, the data recorded represent the sky convolved 
with a point spread function (PSF). The PSF shape de- 
pends not only on the mirror and aperture (window) of 
the telescope, but also often on the atmosphere. Even 
in an ideal instrumental setup, images are blurred due 
to secondary peaks of Rayleigh diffraction. In modern 
astronomy, spatial resolution plays a pivotal role in deci- 
phering morphology and obtaining accurate photometry 
a nd henc e it is important to try to de-blur images. 

IMeinell ([1986D has shown a wide variety of both linear 
and nonlinear restoration algorithms, based on maximum 
likelihood and recursive improvement of images in suc- 
cessive iterations, taking account of the Gaussian and 
Poisson noise processes. Moreover, the additional a pri- 
ori information about the positivity of the true image 
helps to improve resolution. One of the most frequently- 
used methods is the above-mentioned L-R deconvolution 
scheme, which conserves flux. A key property of this 
algorithm is that it converges toward the most likely so- 
lution of the PSF and image intensity distribution. In 
the literature, there exist several other techniques for 
image restoration. A detailed description of the Maxi- 
m um Entropy Method (MEM) of deco nvolution is given 
in iNaravan fc Nityanandal (|1986f l and iBrvan &; Skillinel 
(1980). A Maximum Correlation Method (MCM), which 
maximizes the correlation between adjacent pixels, and 
simultaneousl y improves ima g e reso lution has been been 
prescribed by lAumann eTall IpQCh for improving IRAS 
survey maps. This algorithm is a lso an e xtensi on of 
the L-R solution as described in iMeineli (|1986f l and 
iVelusamv et al.l (|2008f ). The latter authors adopted the 
MCM technique for deconvolution of Spitzer images, im- 
proving the resolution by a factor of about three. 

These techniques generally require that the PSF be 
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previously determined to high accuracy. When the PSF 
is only poorly known, an iterative blind deconvolution 
technique can be used for image restoration, where both 
the data and the PS F are estimated simultaneously in 
the s uccessive steps (jFish et al.l 1995: T sumurava et al.1 
fl99l . 

2.1. Two Deconvolution Algorithms 

In the classical approach to the imaging process, the 
data image D is expressed by 

D{x,y) = I{x,y)*P{x,y)+N(x,y), (1) 

where P{x,y) is the PSF of the telescope, I(x,y) is 
the (unknown) unblurred sky image, the symbol * de- 
notes the convolution operator, and N{x,y) is the noise 
added to the data while imaging. Recursive formulae 
can be obtained from the maximum likelihood solution 
for two different classes of noise distributio n, namely the 
Poisson and the Gaussian noise processes (|Meinell [19861 : 
iVarosi fe Landsmanl |1993f h From equations (31) and 
(64) of iMeinel ([19861 ) we write dire ctly the solution for 
the P oisson noise case in the form (|Varosi fc Landsmanl 
1993): 



/new *o. 



Id 



D 



old 



*p 



*p 



/new = /old +[[D- (/old * P)\ * P 



and that for Gaussian noise as: 



(2) 



(3) 



Here P is the reflected PSF, P(x) = P(-x). For q = 1, 
equation ([2]) reduces to the L-R algor ithm obtained in- 
dependent ly from Bayesian statistics (jRichardsonl 1 1 972t 
[L^icvl[l97l . The Poisson noise process, by its nature, ex- 
cludes any negative solutions whereas the Gaussian noise 
process allows negative solutions, subject to the choice 
of q in equation ([3]). 

Ideally, after an infinite number of iterations, the out- 
put map should converge to the maximum-likelihood so- 
lution. However, in practice, after a finite number of 
iterations in which a close to maximum-likelihood solu- 
tion has been achieved, the smoothn ess of the map starts 
deteriorating (jVelusamv et al.l 120 08). There are no def- 
inite and gene ra lized stoppin g criteria for avoiding this 
([Prasad! [20021 ). iLucvl (|1974h recommended a stopping 
criterion based o n a g oodness-of-fit test. However, sub- 
sequently ILucvl (]1992f ) noted that a higher number of 
iterations might be needed for images with a large num- 
ber of pixels, and he discouraged the use of the former 
method. 

2.2. Implementation on BLAST maps 

BLAST05 raster-scanned targeted regions in the 
Galactic Plane, having different areas ranging from 
3 deg 2 to 10 deg 2 . From the time-ordered data, maps 
were produced using th e optimal map-maker SANEPIC 
([Patanchon et all [2008). The combination of high scan 
speed and low 1// knee, cross-linking where available, 
together with common- mode removal in SANEPIC, pro- 
duces maps retaining diffuse low spatial frequency emis- 
sion. However, preprocessing of the time-ordered data 
to remove very low frequency drifts, plus using a low 
pass cutoff, makes the SANEPIC map average zero - the 



DC level is not known. The iterative L-R deconvolu- 
tion algorithm for Poisson noise (eq. [2]) converges only 
for true positive images. In order to satisfy the positiv- 
ity of the initial data image D(x,y), we add a constant 
level to the whole map (a deconvolution operation on a 
flat map preserves the initial map in its output). The 
value of the constant is estimated by pixel-pixel correla- 
tion of the BLAST images (smoothed to 4') with the cor- 
res ponding 100 p,m IRIS image, IRAS data reprocessed 
by iMiville-Deschenes fc Lagachd (|2005[ ). We carried out 
deconvolutions with different values of the constant, find- 
ing the solution to be robust. After the deconvolution 
operation we subtracted the constant from the resulting 
images. 

Special care has been taken to deal with the edges of 
the scanned region, padding the external area at a value 
equal to the average around the edge of the map. The 
convolution operations in the iterative updates were car- 
ried out by fast Fourier transform (FFT), and so to main- 
tain a periodic boundary condition and smoothness we 
apodized the entire map edge with a sine function. 

The maps can be produced at various pixelizations, 
usually 15", but also 10" or 20". This is sufficient to 
sample the degraded PSF and to sample and recover 
diffraction-limited information even in the 250 /jm chan- 
nel. Our experiments showed that the deconvolution is 
robust to different pixelizations in this range. Note that 
the corrupted PSF is not azimuthally symmetric. There- 
fore, at each passband a synthetic PSF appropriate to 
the particular scan pattern and coverage of a field is 
constructed from the well-samp led (10" grid) telesc ope- 
frame PSFs of each bolometer (|Chapin e t al. 2008]). In 
order to avoid artefacts at large angular scales we have 
apodized the boundaries of the PSFs with a cos 2 taper 
from a radius of 4f5 to 5'7. 

For the actual processing we have used the 
IDL-implemented ' Max likelihood.pro' routine by 
IVarosi fc Landsmanl (|l993). which includes both meth- 
ods of iterative update, appropriate to Gaussian or 
Poisson noise. Probably neither is an accurate noise 
model for the BLAST Galactic maps, which incorporate 
a fluctuating cirrus background, bright cirrus structures, 
and strong compact sources. Empirically, the largest 
errors in the maps are near the bright sources, this 
arising from small pointing errors, the asymmetric 
PSF, and the particulars of the map-making process 
itself. Furthermore, in the reconstruction step implicit 
in the deconvolution (eq. [1]), whenever the synthetic 
PSF is not a perfect representation there are larger 
errors in predicting D where / is brighter. Thus of the 
two, the Poisson algorithm is expected to work better 
because it gives less weight where there are strong point 
sources, i.e., where there are larger errors in the maps. 
Although the noise is not actually Poisson and the 
solution might not be optimal, it is always consistent 
with the data through equation (flj. Note that a similar 
situation arises in the HIRES processing of IRAS data, 
where the noise mod el is likely a complicated hybrid 
(|Aumann et al.lll990D . We have performed experiments 
with both noise models and, as anticipated, the Poisson 
alternative was found to perform better in enhancing 
the structure in the image. 

As shown in Figure [1] L-R deconvolution progressively 
enhances an image, but at the expense of background 
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Fig. 1. — Example of BLAST image deconvolution by the L- 
R algorithm (Poisson noise), at 250 fim for the IC 5146 region. 
Upper left: unprocessed image where compact sources have the 
characteristic imprint of the corrupted BLAST PSF. Upper right: 
deconvolved image after 8 iterations. Lower left and lower right: 
after 64 and 256 interations, respectively. The typical size of the 
compact sources in the final map is about 1'. 

noise amplification at high frequencies after a large num- 
ber of iterations. This noise can be suppressed by con- 
volving the deconvolved map with a restoring beam G 
of FWHM close to the diffraction limit. Mathematically 
this step of restoring the image might be achieved in two 
distinct ways: 1) convolve each side of equation ([T]) with 
G, and treat / * G as the unknown sky map to be found 
by the iterative solution of this new equation; or alterna- 
tively 2) as the effective PSF use (P/G)(x,y), where X 
denotes the Fourier transform of X. We have compared 
the results for these two alternatives with the results of 
the standard deconvolution subsequently smoothed by 
the restoring beam and have found that L-R deconvolu- 
tion gives robust results for all three cases. The images 
shown in this paper have not been restored. 

The reason for this is that despite the complicated 
and asymmetric PSF and the noise in the maps, ap- 
plication of L-R deconvolution removed the worst ef- 
fects of the corrupted beam, thus improving significantly 
the definition of both compact sources and diffuse struc- 
tures. Compared to the results of the direct inversion, 
the background obtained is much smoother, although 
ringing around bright sources has not been completely 
removed. Such ringing is not seen in Figure [1] where 
there are no very bright sources, but can be se en in the 
ma ps of Aquila (iRivera-ingraham et al.ll2010D . Cygnus 
X (|Rov et al.M2010af l. and K3-50 (below), each of which 
contain significantly brighter sources. This same effect is 
shown in the simulations below. 

This ringing is the principal artefact in the deconvolved 
images, as was also found to be the case in using the 
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HIRES MCM algorithm to produce the Infrared Galaxy 
Atlas (IGA) fro m IRAS scans of the G alactic plane 
(|Cao et al.Ml997D . iVelusamv et~aTl (|2008t ) reduced the 
ringing artefacts in deconvolutions of Spitzer data by first 
subtracting the background around the targeted com- 
pact sources. However, the large maps from the BLAST 
surveys are more complex, having large-scale structures 
of different brightness on which stronger sources are su- 
perimposed. Therefore, it is not practical to reduce the 
background everywhere to be close to zero. 

These artefacts are on the scale of the corrupted PSF, 
and consist initially of a ring-shaped depression with 
an outer enhanced ridge. In successive iterative steps 
the artefact evolves, developing a finer series of depres- 
sions and ridges of smaller amplitude spreading from the 
source out into the rest of the image. The average bright- 
ness enclosed within the ringing pattern is very close 
to the local sky background. The intensity ratio of the 
brightest, smallest ring to the deconvolved central peak 
depends somewhat on the local background, but is typi- 
cally no larger than 0.1%. Nevertheless, near very strong 
sources it is quite obvious, adding uncertainty in the es- 
timation of the sky background level near bright (often 
crowded) sources and limiting the detectability of nearby 
faint sources. 

3. TESTS AND GROUND TRUTH 
3.1. Simulations 

We have performed simulations to check the robust- 
ness of the L-R deconvolution scheme. Fake sources were 
constructed by convolving normalized narrow Gaussians 
with the BLAST PSF and then multiplying them with 
random flux densities. Gaussians with intrinsic FWHMs 
ranging from 0.5' to 2' were tested. These simulated 
sources were inserted at different locations in the map 
such that the sources were well separated from each 
other. During the BLAST map-making process, a vari- 
ance map is obtained from a combination of the noise 
in the time stream data and the map coverage by the 
bolometers in the array. We used this variance estimate 
to add noise to the simulated map, but in these tests we 
did not take into account the background structures and 
fluctuations from the diffuse cirrus. 

We then deconvolved the simulated map using the PSF 
and observed the convergence of the measured FWHMs 
of the sources as a function of iteration, as shown in Fig- 
ure [21 In these tests we were able to recover the intrinsic 
size of the sources, convergence occurring after about 64 
iterations. In Figure [3] we have plotted the average frac- 
tional flux densities recovered at different iterations of 
the L-R deconvolution. For isolated sources with various 
FWHMs about 98% of the input flux density is recov- 
ered. It is possible that the remaining flux density is lost 
due to the ringing artefact near the source, for which a 
Gaussian profile used in the fitting is not an accurate 
model. 

Figure 0] shows the performance of L-R deconvolution 
in a simulation in which two compact sources with rela- 
tive flux density 10 are placed only 1 '5 apart, well within 
the corrupted PSF. In the blurred map it is hard to dis- 
cern a faint second source hidden within the brighter 
source. However, L-R deconvolution makes it clear. Fit- 
ting a double Gaussian to the deconvolved map, we have 
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Fig. 2. — Measured FWHM of simulated sources at successive 
iterations of L-R deconvolution. Filled circles, triangles, open cir- 
cles, and crosses represent sources of intrinsic size Of 5, l'.O, 1.'5, 
and 2.'0, respectively. 

retrieved the initial flux densities and sizes within 96% 
and 99%, respectively. The positions of the recovered 
sources are accurate within ~ 1". 

3.2. L-R Performance on Compact Sources 

Application of the L-R scheme to BLAST05 maps of 
the Aquila and Cygnus X regions has already shown im- 
pressive improve ment in resolving confused sources in 
crowd ed regions ivera- Ingraham ct al. 2 010tlRoy et al.l 
l2010a|) . thereby enabling the preparation of deeper cat- 
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Fig. 3. — Measured fractional flux density averaged for compact 
sources with various FWHMs, at successive iterations of L-R de- 
convolution. 
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Fig. 4. — Left panel: two simulated sources of relative flux density 
10 placed 1'5 apart. Right panel: deconvolved map showing the 
two resolved sources. The brightness scale has been chosen to 
highlight the faint ring artefact. 

alogs. From surveys at both longer and shorter wave- 
lengths there was abundant ground truth for the fainter 
sources recovered. To demonstrate the L-R performance 
on a field of compact sources we have selected the region 
near K3-50, which is a young star- forming site contain- 
ing a g roup of H II regi o ns. We make use of images from 
MSX (|Mill et all 11994 lEgan et all 11998) at 8 fj,m, the 
IGA at 60 and 100 um, and the C GPS 21-cm radio con- 
tinuum survey (jTavlor et al.ll2003j ). We also make use of 
spectral line cubes of 12 CO from the Five College Radio 
Astron omy Observatory (F CRAO and of H I from the 
CGPS (jTavlor et, alll2003h . The science results for K3-50 
are deferred to § HI 

Figure [5] shows details of the K3-50 region. The up- 
per left panel is the original SANEPIC 250 ^m BLAST 
map, clearly showing imprints of overlapping corrupted 
PSFs. The upper right panel shows how the L-R decon- 
volution operation has improved the resolution, enabling 
the detection of faint sources which were otherwise hid- 
den. The compact sources have an angular resolution of 
~ 40", which corresponds to a spatial extent of about 
1.7 pc at a distance of 8.7 kpc (jPeeters et al.ll2002l h The 
overlaid tria ngles identify the r adio source components 
obtained by ISamal et all (|2010f ) from 1280 MHz obser- 
vations, showing that there are multiple sources within 
the two bright BLAST sources, still unresolved at the 
improved BLAST resolution. Their names are indicated 
in the lower right panel, which also shows radio emission 
and CO contours. The lower left panel shows the same 
region imaged by MSX (8 fixn), also at relatively high 
resolution. Again, there is good correspondence with 
BLAST. 

K3, called IRASB by ISamal et ail (|2010l) . coincides 
with the point source IRAS 19597+3327A. At its posi- 
tion there is a CO peak at —23 km s^ 1 , the same velocity 
as for the much brighter K4 and K5 sources, and in the 
CGPS H I emission there is an absorption reversal at 
the same velocity. This indicates that K3 is related to 
this complex and is not a galaxy, as it has been alterna- 
tively classified (2MASSX J20013735+3335282). K2 is 
an asymmetric ring at 8 /im and has no radio counter- 
part. Kl is on an arc of emission seen in the BLAST 
and MSX 8 ^tm maps, apparently related to the bubble 
of extended H II emission and its confinement. How- 
ever, there is no prominent counterpart in either MSX 
or IRAS images. There is no 21-cm radio continuum 
counterpart either, all indicating that the source is cold 
and unevolved (see below) . There is no confirming infor- 

18 http:/ /www. astro. umass.edu/~fcrao/ 
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Fig. 5. — Upper left: raw BLAST05 map of K3-50 region at 250 fim, showing that crowded sources remain unresolved. Upper right: same 
region after L-R deconvolution. BLAST compact sources are shown by the circles. The recovered FWHM of K4 (K3-50A) is 36". Triangles 
show the positions of radio sources from[Samal et al. (2010), with labels A to F (K3-50A to K3-50F) in the lower right panel. Lower left: 
8 /im MSX image. The arc-like feature is also seen in the BLAST image. Lower right: solid contours from 21-cm radio continuum (CGPS) 
overlaid on the deconvolved 250 fim BLAST image. The dash-dot contours are the 12 CO line emission (FCRAO) integrated over —30 to 
— 16 km s _1 . Note that the CO map coverage does not extend to Kl and K2. 

mation from CO, since unfortunately, Kl (and also K2) is 
outside the area covered by that emis sion-line survey; the 
low resolution integrated CO map of iDame et al.l (2001) 
hints at an extension in this direction. 



3.3. L-R Performance on Diffuse Structures 

In this section we discuss the effectiveness of L-R de- 
convolution for recovering the morphology of extended 
structures. BLAST05 mapp ed a diffuse field towa rd the 
Cas A supernova remnant (|Sibthorpe et al.lfeOlOP ). The 
map was scanned in only one principal direction, but 
nevertheless the L-R deconvolution has restored numer- 
ous elongated diffuse structures. We have noted how the 
contrast of these structures increases with successive it- 
erations. However, ultimately the background noise be- 
comes amplified into non-physical small-scale structure, 
and so intermediate iterations produce the best com- 
promise for studying the details of the diffuse emission. 
From the appearance of the maps, it was judged that 32 
iterations is optimal. Subsequently, this field has been 



imaged by HerscheJ^ (Barlo w et al.ll2010h . providing the 
desired ground truth to verify the efficacy of the L-R de- 
convolution. For comparison, the 250 /im SPIRE image 
was smoothed with a 40" Gaussian, and regridded to the 
20" pixelization of the BLAST image. Figure [5] shows 
this SPIRE image overlaid with contours from the cor- 
responding deconvolved 250 fim BLAST map after 32 
iterations. We have again examined the higher-iteration 
images and confirmed from this direct comparison that 
for the signal to noise of this BLAST survey, 32 iterations 
is about optimal. 

There is also good evidence for recovering diffuse struc- 
tures in the fields being studied in this paper. The MSX 
8 jum image of the K3-50 region in Figure [5] shows emis- 
sion (probably from PAHs) along an extended shell-like 
structure, labelled "arc." Related structure is also seen in 
the radio emission, suggesting a bubble and surrounding 
PDR The deconvolved BLAST05 map reveals that the 

19 http:/ /herschel.esac.esa.int/Science_Archive.shtml 
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Fig. 6. — SPIRE 250 /an image of the Cas A field convolved 
with a 40" Gaussian and reprojected on the 20" BLAST grid. 
The contour overlay is from the corresponding BLAST deconvolved 
map, after 32 iterations. After dcconvolution, the BLAST image 
faithfully reveals the morphology of the diffuse structures. 

same structure is also being traced by dust continuum 
emission, though the effect of the map edge on the lower 
right is apparent too. Furthermore, there is a faint ridge 
of submillimeter dust emission extending to the lower 
left of source K5 (K3-50C) which is traced by the CO 
contour. 

Within the IC 5146 held is an elongated filamen- 
tary molecular cloud structure commonly known as the 
"Northern Streamer." Figure [7] shows a S CUBA 850 urn 
dust map of a targeted section in IC 5146 (|Kramer et al.l 
l2003t iDi Francesco et all I2008D , overlaid with contours 
from the 500 fim deconvolved image. This shows that 
the L-R operation has preserved both large-scale struc- 
tures and smaller-scale fragments along the ridge. A dis- 
cussion of the astrophysics for this region is deferred to 

m 

4. SCIENCE RESULTS FOR K3-50 

K3-50 is a well-known group of compact H II regions, 
namely K3-50A-F, within the star-forming complex W58 
(|Israellll976t [Howard et al.l[l996t ISamal"et^l20Tol) . The 
warm dust surrounding the most prominent ionized re- 
gion is bright in both far-infrared and submillimeter 
emission; the correlation with the CGPS 21-cm radio 
continuum emission is seen in Figure [5] 

The radio sources K3-50A and K3-50B are well cen- 
tered inside the BLAST clump K4 ( spatial exten t ~ 1.7 
pc), whereas K3-50C is within K5. iSamal et ail (|2010f ) 
further resolved K3-50C into two sources (CI and C2). 
Moreover, this young star-forming region contains em- 
bedd ed infrared star clusters, namely K3-50D and K3- 
50B (jBica et al.ll2"00l . 

In the overlap region between the BLAST survey and 
the FCRAO 12 CO molecular line survey there is good 
correspondence between peaks of submillimeter and CO 
emission, which provides further insight into global mor- 




93.7 

Galactic Longitude 



Fig. 7. — Contours from the deconvolved BLAST 500 (im image 
overlaid on the SCUBA 850 fim map of the "Northern Streamer" 
in the IC 5146 region. BLAST contours trace small-scale sub- 
structures in the SCUBA map along the filament, again showing 
the ability of L-R deconvolution to restore morphological detail 
within diffuse emission. 

phology and gas dynamics in the region. The CO con- 
tours on the lower right panel of Figure [S] show that K3, 
K4, and K5 are part of the same cloud (also s ee § 13. 2[) . 

In our preliminary work ([Truch et al.l [20081 ) we mea- 
sured flux densities for K3-50A and K3-50C (K4 and K5) 
by simultaneously fitting a model of two 2-dimensional 
Gaussians convolved with the corrupted BLAST PSF, 
where sizes, positions, and amplitudes were free param- 
eters of the fit. In our present analysis, with the well- 
resolved deconvolved map, we fit multiple Gaussians with 
a linear background model to extract flux densities for in- 
dividual sources. Table Q] summarizes the flux densities. 

In order to achieve a broader multi- wavelength descrip- 
tion of the SEDs, we have extracted flux densities at 60 
and 100 /jm from the IGA maps. For longer wavelength 
data poin ts we have used the inte grated flux densities re- 
ported bv lThompson et al.l (|2006f l at 450 and 850 /im for 
K4. Due to the restricted coverage of the SCUBA imag- 
ing, flux densities at 450 and 850 /im are not available for 
the other sources. We have also used MSX photometric 
data in deriving the bolometric luminosity; however, we 
have not used these data in the simple fit to the long- 
wavelength SED. 

The SED-fitting model is a simplified single- 
temperature modified blackbody with a fixed dust emis- 
sivity index of /3 = 1.5. To obtain th e mass, we have used 
equation (2) of IRov et al.l (|2010af) . adopting the same 
250- [im dust opacity kq = 10 cm 2 g _1 and dust to gas ra- 
tio r = 0.01. We use the given distance of 8.7 kpc. From 
the H I absorption velocities and the CO emission, the 
sources are in a weak "Perseus arm" feature (in I — v or 
b — v diagrams) beyond the tangent point (at ~ 5.5 kpc), 
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TABLE 1 

SED Best-Fit Parameters of BLAST Sources in the K3-50 Field 



ID 


BLAST Source 


5250 


£350 


S500 


SlOu 


5*60 


T M L L hol 




Name 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(K) (10 2 M ) (10 4 L ) (10 4 L ) 



Kl 


J200054 - 


h 333055 


105 ± 10 


61 ±5 


23 ±3 


70 a 


14 a 


20 ± 1 


28 ±5 


1 ± 1 




K2 


J200112 - 


h 333802 


96 ±9 


46 ±3 


13 ±1 


277 


175 


33 ± 1 


6± 1 


3±1 




K3 b 


J200137 - 


h 333454 


144 ± 20 


74 ±8 


20 ±8 




308 


34 ± 1 


10 ±2 


5±1 




K4 


J200146 - 


h 333214 


1468 ± 19 


576 ±9 


245 ±6 


11600 


12000 


44 ±4 


68 ± 14 


152 ± 54 


192 


K5 


J200154 - 


r 333343 


809 ± 25 


357 ± 10 


184 ±6 


3300 


3100 


35 ±2 


55 ± 14 


36 ±4 


39 



a Data used as upper limit in SED fit 
b IRAS 19597+3327A 



T 




Wavelength [/im] 



Fig. 8.— SED for K3-50A (our K4). Flux density multiplied by 
v/v250 plotted along the y axis. The solid line is the best fit to the 
modified blackbody with a fixed dust emissivity parameter of f3 = 
1.5. The grey lines represent an envelope of 68% confidence level 
obtained from Monte-Carlo simulations. The best-fit parameters 
are given in Table \T\ Data at wavelengths lower than 60 fira are 
not included in the fit but are used in estimating Lbol- 

but not beyond the stronger emission at —60 kms . 
Table [1] summarizes the physical properties, namely tem- 
perature (T), luminosity (£), and mass (M) for the indi- 
vidual BLAST sources. The reported 1-er uncertainties of 
the physical parameters ar e obtained from Mon te-Carlo 
simulations as described in IChapin et al.l (|2008f) . 

Figure [H] shows an example of an SED fit, for K4. A 
best-fit temperature of 43 ± 4 K and total mass of (7± 
I)xl0 3 M(T) are obtain ed from the parameters of the fit. 
IQkamoto et al.1 (|2003l ) have found a star cluster associ- 
ated with K3-50A (the dominant rad io source inside K4 , 
also an ultra-compact H II region, IKurtz et al.lll99l . 
Our derived luminosity is (1.5 ± 0.5) xlO 6 L Q . This 
clump lu minosity is equ ivalent to a ZAMS sp ectral type 
of ~Q 4 (TPanagial 119731) . which agrees with iSamal et al.l 
(2010), who assigned a spectral type earlier than 05, 
based on near-infrared measurements. However the high 
luminosity supports the possibility of multiple ionizing 
stars in the stellar cluster. In fact, for K4 it is impor- 
tant to note that we measure dust emission whose energy 
comes collectively from K3-50A and K3-50B. 

The position of a source on the L-M p lane provides a 
rough estimate of its evolutionary stage (Molinari et al. 
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Fig. 9. — BLAST compact sources in the K3-50 region in the L-M 
plane. Solid diagonal lines are loci of constant T or L/M. Dot-dash 
curves orthogonal to these are for constant 250 fim flux density. 
Dash and dot-dash lines denote sources powered by accretion and 
nuclear burning, resp ectively, as derived empirically in Fig. 9 of 
IMolinari et~all 1200ST ). 

2008). The relationship of the evolutionary sequence to 
the underlying energe tics powering the d ust emission is 
discussed in detail by IRov et all (|2010aj ) . Figure |9] also 
confirms that K4 contains one or more zero age main 
sequence stars, deriving luminosity mainly from nuclear 
burning, and hot enough to cause copious ionization. 

In a similar fashion, for K3- 50C (K5) we obtain ed a 
temperature 35 ± 2 K, which iTruch et all (j2008f ) un- 
derestimated due to confusion in the maps combined 
with the uncertainties in the assumed model. The en- 
ergetics of K5 are primarily influenced by the two radio 
sources inside the clump, namely K3-50C1 and K3-50C2 
(jSamal et a l. 2010). Compared to K3-50A, the dust tem- 
perature is somewhat lower, suggesting that K3-50C is 
deeply embedded inside the molecular cloud, is less ener- 
getic, and is perhaps somewhat less evolved. Integrating 
the SED for K3-50C, we obtain a luminosity of (3.6 ± 
0.4) xlO 5 L w hich is equivalent to -06 on the ZAMS 
(|Panagialll973D . The clump mass is (5.5 ± 1.5) xlO 3 M Q , 
placing it in the appropriate position in the L-M dia- 
gram. 

K3 is located between the two radio sources K3- 
50E and K3-50F, to its West and East, respectively 
(see Fig. [5 |) It s luminosity is (5 ± l)xl0 4 L Q ; 
ISamal et all ((2010) underestimated the luminosity (~ 
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2.5 xlO 3 L©) due to the absence of FIR/submillimeter 
coverage. The luminosity and mass of (10 ± 2) x 10 2 M Q 
position K3 in the L-M plane as a "class I" object, 
though with possibly some power coming from accretion. 
Its equivalent single-star ZAMS spectral type is ^08.5, 
indicating the possibility of ionizing its surroundings. No 
significant radio continuum emission peak is detected ei- 
ther at 1420 MHz (see Fig. El) o r above the 5 mJy level 
at 1280 MHz (jSamal et alJ l2010). Possibly there is radio 
self-absorption, or the ionizing radiation is absorbed by 
dust in a dense envelope. Alternatively, the luminosity 
may come from many slightly lower mass stars, which, 
being cooler, would produce collectively less ionizing ra- 
diation. 

A general introduction of Kl and K2 was already given 
in § 13.21 An accurate distance estimate for these sources 
is not available. For our present analysis we have as- 
sumed a distance of 8.7 kpc, similar to the K3-50 region. 
However, note that the assessment of the evolutionary 
stages from the L-M diagram is not affected by the dis- 
tance uncertainty; the position on the L-M plane simply 
shifts diagonally along a line of constant LjM or tem- 
perature. The K2 clump is associated with emission in 
the IRAS 60 and 100 /mi bands and has extended MSX 
8 /im emission as well. It is relatively hot, with a tem- 
perature of 33 ± 1 K. The L-M plot suggests it is a less 
energetic version of K3, though the MSX morphology is 
distinctly different. Kl, projected at the outskirts of the 
K3-50 region, with no YSO or radio counterpart, is the 
coldest BLAST clump in this region, at 19 ± 1 K. In 
the L-M diagram it is at an earlier stage of evolution, 
beginning to be powered by accretion. 

5. SCIENCE RESULTS FOR IC 5146 

The IC 5146 molecular cloud comple x in Cygnus has 
been widely studied i n the optical, IR flHerbig fc Dahml 
2002) submillimeter (Kramer et alj [2003) and molecular 
lines (Dobashi et al.lll992f). Using infra red color excesses 
in the JHK bands. lLada et all (| 19941 ) produced an ex- 
tinction map, i.e., the spatial distribution for the dust 
column density. They also surveye d in the molecular 
line emission of 13 CO, C 18 0, and CS. lHerbig k. ReipurthI 
(2008| provide a summary of the o bservational progres s 
achieved in this star-forming region. lHarvev et atl ()2008l ) 
have surveyed the IC 5146 region with Spitzer to study 
the properties of young stellar objects. Distance esti- 
mate s range from 460 pc to 1.4 kpc (see lHarvev et al.1 
120081 and references there in). We have adopted 1 kpc 
from lDobashi et aLl (11992]), which is close to the 950 pc 
used bv lHarvev et al.l ((2008) . Proximity is an advantage 
for probing the workings of star-forming regions. At this 
distance, BLAST can resolve spatial structures of about 
0.3 pc. Our BLAST05 observation (Fig. [10]) targeted the 
filamentary structure of the "Northern Streamer," one 
of the densest molecu lar clouds in IC 5146 ("cloud C" 
of lDobashi et al.lll992f ) , and reportedly the most massive 
(^2.2 xlO 3 M Q ). The field also contains the dark cloud 
L1030 and L1031 to the South- West. 

5.1. Compact Sources 

Figure [TO] shows the locations of ten BLAST compact 
sources in the IC 5146 region. Six of these, with IRAS 
counterparts, were identified as protostellar candidates 
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Fig. 10.— BLAST map of IC 5146 at 250 fj,m, showing relation- 
ships to other cataloged objects. Circles of radius 1' locate the 
BLAST compact sources. Large squares correspond to the peak 
positions of the ve locity-integrated 13 C O (J = 1 — > 0) molecular 
line emission map (Dobashi ct al. 1992). YSOs from Harvey et al. 
(2008) are shown by small circles, stars, triangles, and crosses, to 
distinguish Class I, flat, II, and III spectra, respectively; the plus 
symbols represent YSO candidates for which there is no classifica- 
tion. The dust temperature of the c old fi lament within the outlined 
do t-dashed box i s 11.7 ± 0.4 K (§ [572}. Four cores characterized 
bv I Kra mer ct al. (2003;) are marked with large crosses. Some ad- 
ditional cold and starless clumps identified by BLAST are marked 
by large diamonds. 

by iDobashi et~al"1 (I1992D. Amon g these, Eliasl-12 (our 
IC9) is a T-Tauri star (|Elias 1978 ) , also known t o be asso- 
ciated with molecular outflows (jLevre ault 1983). Follow- 
up surveys with 12 CO (|Dobashi et al.l 119931 I2001D have 
shown that all of these sources have molecular outflows, 
revealing the early stages of star formation in progress in- 
side the embedding clumps seen by BLAST. The squares 
in Figure [TO] correspond to the peak positions of the 
velocity-integrated 13 CO (J = 1 — > 0) molecular line 
emission map. There is no direct correspondence be- 
tween the BLAST dust emission peaks and the 13 CO 
peaks, but the BLAST compact sources are associate d 
w ith the CO clumps ( see Fig. 6 of lDobashi et al.lll992t) . 

lHarvev et all (|2008l ) have detected 60 YSO candidates 
within the region surveyed by BLAST. Using the slopes 
across MIPS and IRAC data they classified 38 objects 
into Class I, flat, II, and III. Thes e relate to progres - 
sively later stages of YSO evolution (Evans et al. 2009). 
Their positions and classifications are given in Figure [TO] 
Along the dense "Northern Streamer" most of the sources 
belong to Class I and flat classifications, i.e., earlier 
stages of evolution. Likewise, there are Spitzer YSO 
counterparts to the ten BLAST sources, these mostly 
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TABLE 2 

Flux Densities of BLAST Sources in the IC 5146 Field 



ID 


BLAST Source 


I 


b 


F250 


^350 


-P500 


F100 


^60 


F 2i 




Name 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 



ICl 


J214558 - 


h 471835 


93.4334 


-4.6663 


97.1 ± 9.0 


75.1 ±7.3 


32.5 ±2.1 








IC2 a 


J214508 - 


h 473305 


93.4842 


-4.3840 


42.5 ± 1.5 


39.7 ± 1.4 


22.4 ±0.8 


29.0 ±0.5 


14.2 ±0.1 


0.7 


IC3 


J214626 - 


h 471757 


93.4844 


-4.7141 


60.1 ± 4.7 


42.3 ±2.3 


24.0 ± 1.0 






0.1 


IC4 b 


J214452 - 


h 474026 


93.5297 


-4.2616 


38.4 ±2.0 


31.5 ±1.5 


19.9 ± 1.2 


24.3 ±0.3 


13.3 ±0.01 


0.8 


IC5 


J214534 - 


h 473559 


93.5703 


-4.3936 


15.0 ± 1.0 


12.1 ±0.6 


8.7 ±0.4 






0.1 


IC6 C 


J214443 - 


h 474635 


93.5769 


-4.1670 


45.1 ± 3.4 


33.5 ± 1.7 


19.8 ± 1.2 


10.1 ±0.3 


2.3 ±0.1 


0.2 


IC7 


J214558 - 


h 473536 


93.6132 


-4.4397 


69.1 ± 6.5 


47.7 ±2.3 


30.0 ± 1.3 


20.0 ±6.2 


1.7 ±0.1 


0.2 


IC8 


J214707 - 


h 473242 


93.7369 


-4.6054 


23.0 ±2.6 


20.6 ± 1.5 


11.4 ± 1.5 








IC9 d 


J214722 - 


h 473221 


93.7659 


-4.6372 


186.6 ±3.4 


121.8 ±2.2 


66.9 ± 1.2 


97.5 ± 10.0 


36. ± 2.2 


2.5 


IC10 


J214759 - 


h 473643 


93.8913 


-4.6470 


45.4 ±2.7 


30.8 ±1.3 


14.7 ±0.8 


12.5 ±1.2 


2.4 ±0.2 


0.3 



a IRAS 21433±4719 
b IRAS 21429±4726 
c IRAS 21428±4732 
d V1735 Cyg, Eliasl-12 

being Class I. These YSOs within the more extensive 
BLAST clumps (0.3 pc) are of relatively low power and 
so the dust is not significantly internally heated and the 
BLAST dust temperatures are all quite low. For exam- 
ple, IC9 contains several YSOs (in the 2MASS if-band 
image there is reflection nebulosity) . At the present sub- 
millimeter sensitivity and resolution, not all YSOs have 
BLAST counterparts. 

We have characterized the physical properties of the 
identified BLAST clumps by fitting SEDs as described 
in § HI Compared to the K3-50 region (dominated by 
warm dust emission), the dust temperature of IC 5146 
is relatively cold, which motivated us to constrain SEDs 
with a spectral index of (3 = 2.0. 

Where there are separable IRAS sources, before per- 
forming photometry on the BLAST sources, we have con- 
volved the respective maps to the 2' native beam res - 
olution of the IGA 100 fj,m images (|Cao et al.l 119971 ) . 
Otherwise we used the better BLAST resolution. Ta- 
ble [2] summarizes the flux densities. All of these sources 
have 250 fim flux densities greater than at 100 /jm, in- 
dicating a relatively cold dust temperature. Constrain- 
ing the SEDs with the longer wavelength data points 
is important, but unfortunately SCUBA has observed 
only a small part of the IC 5146 region (Fig. [7]). This 
includes IC7, IC8 and IC9, for which we obtained in- 
tegrated 850-^im flux densities of 10.6, 1.5, and 5.6 Jy, 
respectively. For the cool sources here we have used only 
the 100 /im flux density in the SED fits, while the 60 and 
24 p,m flux dens i ties p rovid e upper limits, as d iscussed 
in IChapin etaLl (|2008l) and iTruch et al.l (|2008[ ). As an 
example we show the SED for the brightest clump, IC9, 
in Figure ITT1 

Table [3] summarizes the physical quantities T, M, and 
L obtained from the best-fit parameters of the SEDs. In 
this region of IC 5146 the luminosities of the identified 
BLAST sources range from 11 to 170 L Q and the masses 
from 15 to 180 M Q . The temperatures range from 13 to 
16 K, with an average about 14.8 K. 

When a clump is more massive than the Bonner-Ebert 
critical mass, then gravity overpowers the internal (ther- 
mal, magnetic, and turbulent) supporting pressure and 
it undergoes collapse. Our BLAST mass estimates ap- 
pear to be above the critical mass (using the simple 
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Fig. 11.— SED as in Fig.|8] but for IC9 (Eliasl-12) in the IC 5146 
region fit with f} = 2. Arrows for short-wavelength data represent 
3-<T upper limits above the measured flux density. The best-fit 
temperature for the cold envelope is 16 K. 



TABLE 3 

SED Best-Fit Parameters of BLAST 
Sources in the IC 5146 Field 



ID 


T 


Mtot 


L 




(K) 


(M ) 


(L©) 


ICl 


12.7 ± 1.2 


171 ± 50 


48 ±16 


IC2 


16.3 ±0.6 


35 ±8 


43 ±4 


IC3 


14.9 ±0.1 


64 ±6 


47 ±5 


IC4 


16.1 ±0.7 


33 ±9 


38 ±4 


IC5 


14.9 ±0.1 


15 ±2 


11 ± 2 


IC6 


13.6 ±0.3 


73 ± 11 


31 ± 2 


IC7 


14.0 ±0.3 


106 ± 10 


53 ± 3 


ICS 


15.0 ±0.1 


25 ±3 


19 ± 2 


IC9 


15.7 ±0.3 


176 ± 18 


171 ± 12 


IC10 


14.3 ±0.3 


55 ±7 


31 ± 2 
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prescription described in lKerton et al.ll2~001l ). suggesting 
that the compact sources are gravitationally unstable. 
Indeed, there are already some Y SOs that are generat- 
ing outflows. lDobashi~e t al. (200T|) showed how physical 
processes relating to mass, momentum, and energy car- 
ried away by the outflows have direct consequences for 
the stability of the parent cloud and as well influence 
the evolutionary dynamics of the embedded stars. These 
outflow parameters have a correspondence with the bolo- 
metric luminosity of the accreting source (Dobashi et al. 
120011 ). Certainly feedback pr ocesses have al r eady begun 
in this molecular cloud, and Dobas hi et al.l (|1992[ ) have 
suggested that outflows have played an important role in 
supporting the parent cloud from collapsing. 

5.2. Cold Central Filament 

Apart from the above-mentioned compact sources, 
BLAST observes abundant substructure in the North- 
ern Streamer reg ion. The infrared extinction map by 
lLada et "all (|1994f) has also revealed a distribution of high 
column density clumps throughout this filament. There 
is a good correlation of Ay wit h the 12 CO an d 13 CO 
molecular line emission maps bv iDobashi et all |l992). 
Not unexpectedly, our dust emission map also has a strik- 
ingly tight correlation with Ay (and CO which might not 
be a perfect tracer of column density at large spatial den- 
sity) , and with 850 /im emission, a s shown in Fig ure [7] 

With earlier SCUBA observations lKramer et ail ([2003) 
studied a narrow 14' by 2.'5 region extending from IC7 
(only at the boundary of their map) to IC8 (included). 
From the SCUBA colors they found dust temperatures 
ranging up to 20 K at the outskirts and between col- 
umn density peaks and down to 10 K in "cores" (coin- 
cident with peaks in Ay), embedded condensations that 
effectively shield the interstellar radiation from which the 
dust derives its power in the absence of a YSO. Note that 
along the central filament studied there are no embedded 
YSOs (Fig. [TU]), except at the two ends. 

We have estimated a global dust temperature for the 
portion of the streamer within the box shown in Fig- 
ure [TUl by fitting a relative SED to the pixel-pixel corre- 
lation slopes of the image data with respect to 250 /im. 
We have used the 850 ^m SCUBA map to constrain the 
spectrum at long wavelengths. The filament is not promi- 
nent in IRAS bands because the dust is cold. Figure [T2l 
shows the SED. With a dust emissivity index /3 of 2.0 the 
best- fit temperature is 11 .7 ± 0.4 K, which agrees well 
with Kramer et all (|2003| ) . 

5.3. Cold Starless Clumps 

Because the BLAST bands are sensitive to colder 
dust emission we can detect cold clumps early in their 
pre-stellar evolution, when their primary ene rgy source 
is the external radiation field ( "stage E" in iRov et all 
l2010al) . Good examples c onfirmed by BLAST are 
the four "cores" identified by Kramer et al.l (|2003l) and 
mar ked by large crosses in Figure[l0] The core numbered 
2 bv lKramer et al.l (|2003f ) is the warmest (18.2 K). This 
is a complex region near IC8 with large temperature gra- 
dients influenced by a YSO visible even to IRAS at short 
wavelengths. However the other three cores are cold and 
apparently starless. 

With the larger spatial coverage of the BLAST sur- 
vey we have identified other cold clumps that have no 
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Fig. 12. — Relative SED obtained from pixel-pixel image cor- 
relation with respect to 250 fi m em ission for the area within the 
dot-dashed box marked in Fig. 1101 With the dust emissivity in- 
dex fixed at 2.0, the best-fit temperature is 11.7 ± 0.4 K. Data 
points represented by circles are not actually used in the fit, while 
the triangle shows the value of unity at the normalizing 250 fim 
point. 

signs of star formation yet. Some prominent examples 
are marked by large diamonds in Figure fTOl at the posi- 
tions of the peaks in the 500 jim. map. Compa ring with 
the visual extinction map of lLada et al.l (|1994[ ). we find 
that these BLAST clumps G93.43-4.28, G93.48-4.30, 
and G93.54— 4.33 have Ay about 10, 20, and 6 mag, re- 
spectively. 

The morphology in this field is complex and variable 
with wavelength, because of the different stages of evolu- 
tion occurring in close proximity. For example, the pro- 
tostar IC4 occurs at the end of an elongated condensation 
that has a cold extension in the direction of G93.54— 4.33. 
Studies of this region will benefit from the improved reso- 
lution of Herschel. With map zero points available from 
Planck, it will be interesting to map the temperature 
struct ure in detail, as has been accomplishe d in other 
fields (|Bernardet al.ll20Tol: IJurela et al.ll20iol ). 

6. CONCLUSION 

The Lucy-Richardson deconvolution algorithm has 
been applied successfully to improve BLAST05 images 
obtained with a corrupted (but known) PSF. This de- 
convolution has enhanced the raw BLAST map of ^3f3 
resolution to near the anticipated diffraction limit. 
This improves the detectability of faint sources; diffuse 
structures are revealed in finer detail as well. We have 
checked via simulations the robustness of applying the 
L-R scheme, especially the aspects of conserving flux, re- 
liability of restoring intrinsic sizes, and performance in 
resolving sources in a crowded field. For the actual de- 
convolved maps of BLAST05 survey fields we have pro- 
vided further ground truth for the improved detail in 
crowded fields and for diffuse structures, by comparing 
with available multi-wavelength high resolution images 
of dust and other tracers. 

We have presented science results enabled by the de- 
convolved maps of two star-forming regions covered by 
BLAST05, namely K3-50 and IC 5146. We were able to 
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resolve three crowded sources in K3-50, namely K3-50A, 
K3-50C, and IRASB, and have also characterized another 
two sources which had previously remained undetected. 

The deconvolved maps of IC 5146 have further shown 
the richness of the field, consisting of both large scale dif- 
fuse structures and compact sources, with a considerable 
range of dust temperature because of the different stages 
of evolution. The compact BLAST sources characterized 
here all have associated YSOs. However, these are not so 
powerful as to heat up the entire dust clump measured 
by BLAST, so that L/M < 1 L Q /M Q and the effective 
dust temperature is rather low; star formation is just 
getting under way. In fact elsewhere there are high col- 
umn density structures that are starless, like within the 



Aumann, H. H., Fowler, J. W., & Melnyk, M. 1990, AJ, 99, 1674 
Barlow, M. J., Krause, O., Swinyard, B. M., Sibthorpe, B., Besel, 
M., Wesson, R., Ivison, R. J., Dunne, L., Gear, W. K., Gomez, 
H. L., Hargrave, P. C, Henning, T., Leeks, S. J., Lim, T. L., 
Olofsson, G., & Polehampton, E. T. 2010, A&A, 518, L138+ 
Bcrnard, J., Paradis, D., Marshall, D. J., Montier, L., Lagache, 
G., Paladini, R., Veneziani, M., Brunt, C. M., Mottram, J. C, 
Martin, P., Ristorcelli, I., Noriega-Crcspo, A., Compiegne, M., 
Flagey, N., Anderson, L. D., Popcscu, C. C, Tuffs, R., Reach, 
W., White, G., Benedetti, M., Calzoletti, L., Digiorgio, A. M., 
Faustini, F., Juvela, M., Joblin, C, Joncas, G., 
Mivilles-Dcschencs, M., Olmi, L., Traficante, A., Piacentini, F., 
Zavagno, A., & Molinari, S. 2010, A&A, 518, L88+ 
Bica, E., Bonatto, C, & Dutra, C. M. 2003, A&A, 405, 991 
Bryan, R. K. & Skilling, J. 1980, MNRAS, 191, 69 
Cao, Y., Terebey, S., Prince, T. A., & Beichman, C. A. 1997, 

ApJS, 111, 387 
Chapin, E. L., Ade, P. A. R., Bock, J. J., Brunt, C, Devlin, 
M. J., Dicker, S., Griffin, M., Gundersen, J. O., Halpern, M., 
Hargrave, P. C, Hughes, D. H., Klein, J., Marsden, G., Martin, 
P. G., Mauskopf, P., Netterfield, C. B., Olmi, L., Pascale, E., 
Patanchon, G., Rex, M., Scott, D., Semisch, C, Truch, 
M. D. P., Tucker, C, Tucker, G. S., Viero, M. P., & Wiebe, 
D. V. 2008, ApJ, 681, 428 
Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 
Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & 

Ledwosinska, E. 2008, ApJS, 175, 277 
Dobashi, K., Onishi, T., Iwata, T., Nagahama, T., Patel, N., 

Snell, R. L., & Fukui, Y. 1993, AJ, 105, 1487 
Dobashi, K., Yonekura, Y., & Hayashi, Y. 2001, PASJ, 53, 811 
Dobashi, K., Yonekura, Y., Mizuno, A., & Fukui, Y. 1992, AJ, 
104, 1525 

Egan, M. P., Shipman, R. F., Price, S. D., Carey, S. J., Clark, 

F. O., & Cohen, M. 1998, ApJ, 494, L199+ 
Elias, J. H. 1978, ApJ, 223, 859 

Evans, N. J., Dunham, M. M., j0rgensen, J. K., Enoch, M. L., 
Merin, B., van Dishoeck, E. F., Alcala, J. M., Myers, P. C, 
Stapelfeldt, K. R., Huard, T. L., Allen, L. E., Harvey, P. M., 
van Kempen, T., Blake, G. A., Koerner, D. W., Mundy, L. G., 
Padgett, D. L., & Sargent, A. I. 2009, ApJS, 181, 321 

Fish, D. A., Brinicombc, A. M., Pike, E. R., & Walker, J. G. 
1995, Journal of the Optical Society of America A, 12, 58 



central filament in the Northern Streamer, for which we 
have obtained an average dust temperature of 11.7 K. We 
have found further examples in the wider BLAST field 
surveyed. 
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